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ABSTRACT 

We collate the evidence for rotation-powered neutron stars that are visible as X-ray 
sources and not as radio pulsars. To date, ten objects have been proposed and one, 
Geminga, has been confirmed as a pulsar by the detection of 4.2 Hz pulsations. Several 
indicators have been used to support the proposition that the X-ray sources are iso- 
lated neutron stars, including high X-ray to optical/radio flux ratios, a constant X-ray 
flux and coincidence with a 7-ray source. Seven of the published neutron star candi- 
dates are located near the centres of supernova remnants, two of them within plerions, 
suggesting that these are young objects (r <20,000 yr). The remaining candidate neu- 
tron stars have no associated supernova remnant and may be older systems, powered 
either by their rotation, like Geminga, or possibly by accretion from the interstellar 
medium. 

Quantitative upper limits exist for the radio fluxes of eight of the ten objects and 
reveal a population at least an order of magnitude less luminous at radio wavelengths 
than known radio pulsars of similar power or age. These could be intrinsically low 
luminosity pulsars, but this implies an overpopulation of neutron stars relative to 
the galactic supernova rate. A simple alternative explanation within the context of 
existing pulsar models is that these objects are pulsars in which the radio beams are 
directed away from Earth. They are still visible as X-ray sources because the weakly 
modulated, surface (thermal) emission, which dominates the soft X-ray emission in 
most young to middle-aged radio pulsars, is radiated in all directions. In the cases 
where hard X-ray or 7-ray fluxes are seen, the beaming explanation implies different 
emission sites for the non-thermal high-energy radiation and the unseen radio beams. 
From the number of candidate neutron stars and radio pulsars younger than 20,000 
years and within 3.5 kpc, the radio beaming fraction of young pulsars is estimated 
to be roughly 50% and certainly much less than 100%. We find the local neutron 
star birth rate to be at least 13 Myr -1 kpc~ 2 . This extrapolates to a galactic rate 
of one neutron star born every ~ 90 years. We conclude that probably all neutron 
stars are born as radio pulsars and that most young, nearby pulsars have already been 
discovered. 
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1 INTRODUCTION 

The Princeton pulsar catalogue (Taylor et al. 1993) now 
contains entries for over 700 radio pulsars, collected from 
many pulsar surveys and targeted searches. In addition to 
the radio pulsars, there is a single entry with no measured 
400 MHz or 1400 MHz flux: the X- and 7-ray pulsar known 
as Geminga (Halpern & Holt 1992, Bertsch et al. 1992). The 
tight upper limits on its radio flux give Geminga a luminosity 
limit several orders of magnitude below those seen in known 
radio pulsars of similar age or power and place this source 
apart from the general population. 

At high energies, there is no obvious difference between 
Geminga and the ~ 20 other pulsars in the Princeton cata- 



logue that have been detected as X-ray sources. Seven, in- 
cluding Geminga, have also been seen by 7-ray telescopes. 
The high energy detections are limited to the very bright- 
est objects and sample a more powerful and nearby set of 
pulsars than the deeper radio searches. In X-rays, the lu- 
minosity is approximately correlated with pulsar spin-down 
power, E, so that E divided by the square of the distance d is 
a convenient measure of the detectability of a radio pulsar as 
an X-ray source. The radio pulsars detected in the ROSAT 
(0.1-2.4 keV) band all have E/d 2 £ 10 34 erg s^kpc' 2 and, 
conversely, nearly all of the radio pulsars fitting this crite- 
rion have been detected. It is important to note that the de- 
tections include not only young, powerful pulsars, but also 
old and millisecond pulsars that are less powerful but are 
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nearby. Taking distance uncertainties into account, the high 
level of X-ray detections therefore implies that most radio 
pulsars are X-ray pulsars. The inverse, that most or all X- 
ray pulsars are radio pulsars, need not be true: Geminga 
is a specific counter-example. Note that the term 'X-ray 
pulsar' is taken here and throughout this paper to mean a 
rotation-powered pulsar observed in X-rays, not a neutron 
star powered by accretion from a stellar companion or from 
a residual accretion disk. 

Geminga is a 350,000-year-old pulsar at a distance of 
only ~ 160 pc. Its proximity has allowed very low lumi- 
nosity limits to be calculated from the limits on its radio 
flux and it has been labelled 'radio-quiet'. However, a low 
flux density does not necessarily mean that the radio lumi- 
nosity is low. It could simply mean that the radio beam is 
not visible from the Earth. In this paper we will define a 
radio-quiet pulsar as a rotation-powered pulsar which has 
not been detected at 400 or 1400 MHz and which therefore 
has a low inferred luminosity at these frequencies. This def- 
inition does not distinguish between explanations based on 
luminosity and beaming, nor does it equate radio- quiet with 
radio- silent. 

Radio-quiet pulsars are extremely difficult to identify 
because photon statistics at high energies make useful pul- 
sation searches impossible for all but the brightest sources. 
The first evidence for a radio-quiet pulsar is usually an ex- 
treme spectrum: bright in X-rays and/or 7-rays but very 
faint at optical wavelengths. While Geminga is the sole con- 
firmed radio-quiet pulsar, nine further objects have been dis- 
covered as unresolved X-ray sources and proposed as radio- 
quiet, rotation-powered neutron stars after fulfilling the first 
criteria for detection of a pulsar. In this paper, we add these 
to Geminga to examine the case for radio-quiet pulsars as 
a class and to distinguish between beaming and luminosity 
explanations for these objects. 

2 THE CANDIDATE PULSARS 

Nine candidate pulsars are listed in Table 1, together with 
Geminga. Each of the candidates is an unresolved X-ray 
source that has not been associated with a compact radio ob- 
ject, leading to the proposal that they are rotation-powered, 
radio-quiet neutron stars. Seven of the objects are within su- 
pernova remnants. Beyond these generalities, the details di- 
verge of the sources and the work done to identify them. The 
table does not include the 'anomalous X-ray pulsars' discov- 
ered in SNRs - 1E2259+586, RXJ1838.4-0301 and 1E1841- 
045 (Gregory & Fahlman 1980, Schwentker 1994, Vasisht & 
Gotthelf 1997) - which are thought to be isolated but pow- 
ered by accretion from a residual disk (van Paradijs et al. 
1995). 

The X-ray fluxes listed in Table 1 provide a useful way 
to estimate the spin-down power E of the candidate pulsars. 
Older radio pulsars detected as point X-ray sources follow 
an approximate trend of L x ~ 3 x W~ 4 E (cf. Seward & 
Wang 1988, Ogelman 1995) which we can invert to estimate 
the spin-down power from the X-ray luminosity. Since the 
trend includes only the detected X-ray pulsars, it is biased 
towards high X-ray luminosities and an E derived from it 
may be underestimated. No allowance for the bias is made 
in this paper because the candidate pulsars are also selected 
as X-ray sources. 



For 3C58 and CTA 1, the plerions provide an inde- 
pendent way to estimate the power of the embedded pul- 
sars. Seward & Wang (1988) found an empirical relation- 
ship between the X-ray luminosity of a plerion and the 
spin-down power of the pulsar powering the plerion. For 
CTA 1, Slane et al. (1997) used this relationship to esti- 
mate a spin-down power of 1.7 x 10 36 ergs _1 . For 3C58, 
similar considerations of the X-ray synchrotron emission 
yield E = (2-4) x 10 36 ergs _1 (Helfand et al. 1995), while 
E — 1.5 x 10 36 ergs _1 is required to explain a radio filament 
near to the candidate pulsar in terms of a shock in the pulsar 
wind (Frail & Moffett 1993). 

The pulsar candidates in CTA 1 and G078. 2+2.1 are 
coincident with 7-ray sources, providing a third way to esti- 
mate the spin-down power. While the relationship between 
7-ray luminosity and spin-down power is still unclear (Fierro 
1995), Brazier et al. (1996, 1997) showed that the 7-ray lu- 
minosities inferred for the two candidates were consistent 
with those of Vela and PSR B1706-44. In both cases, the 
spin-down power derived from the X-ray flux is more than 
an order of magnitude lower than the 7-ray comparison 
would suggest. The estimates from the 7-ray fluxes are listed 
in Table 1. An intermediate spin down power is used in this 
paper for these two objects. 

It would be unwise to treat the candidates in supernova 
remnants as a statistical sample. ROSAT, used to observe 
each of them, has not observed many of the Galaxy's SNRs, 
and identification of a neutron star candidate depends on 
the interests of the observer as well as an object's physi- 
cal properties. In addition, the sensitivity to point sources 
within an SNR varies with the intensity profile of the SNR 
emission, so that objects close to the rim of a SNR may be 
submerged in bright diffuse emission. In RCW 103, for ex- 
ample, the ASCA observation has confirmed the status of 
the point source IE 161358-5055 by separating it spectrally 
from the SNR (Gotthelf et al. 1997). 

In addition to the candidate pulsars within supernova 
remnants, there are a further two candidates, plus the con- 
firmed pulsar Geminga, which are not associated with SNRs. 
The list is short because there is much less to attract the at- 
tention of observers than in the case of an SNR. Without the 
bonus of a supernova remnant, it is also difficult to estimate 
the distances to these pulsar candidates or to guess their 
ages. Nevertheless, the one confirmed radio-quiet pulsar is 
in this group. 

2.1 Candidate pulsars within SNRs 

2.1.1 RXJ0002+6246 (G117.7+0.6) 

A faint X-ray arc was discovered by Hailey & Craig (1995) 
in the periphery of ROSAT PSPC images of the SNR CTB 
1. The arc was proposed to form part of a previously un- 
known shell-type supernova remnant, G117.7+0.6 [footnote: 
The proposed remnant is centred at galactic coordinates 
(l,b) = (117.4,0.3). The position (117.7,0.6) appears to 
be a mis-conversion from celestial to galactic coordinates.] 
at a distance of ~ 3 kpc. An unidentified point source, 
RXJ0002+6246, is clearly visible 10 arcmin from the cen- 
tral position. None of the nearby optical sources was found 
to be a plausible counterpart and there is no radio source 



3 



Name 


F x (0.1 - 2.4keV) 


E 


S400 


Sl400 


SNR 


Common dist 


age 


refs 




(ergem -2 s _1 ) 


(10 30 ergs" 1 ) 


(mjy) 


(mjy) 


name 


name (kpc) 


(10 3 yr) 




RXJ0002+6246 


2 X 10" 13 




? 


? 


G117.7+0.6 


3 


20 


9 


RXJ0007.0+7302 


9 X 10- 14 


7 X 10 6 


< 1.5 


< 0.3 


G119.5+10.2 


CTA1 1.4 


5-10 


1, 2, 16 


RXJ0201. 8+6435 


? 


(2 - 5) x 10 6 


< 2.1 


< 0.15 


G130.7+3.1 


3C58 3.2 


0.8 


3, 7, 16 


IE 0820-4247 


3 x 10" 12 




< 1.7 


< 0.3 


G260.4-3.4 


Pup A 2 


3.7 


4, 14 


IE 1207.4-5209 


2 x 10- 12 




< 1.0 




G296.5+10.0 


PKS 1209-51/52 1.5 


7 


14, 15 


IE 161348-5055 


7 x 10" 13 




< 1.5 


< 0.1 


G332.4-0.4 


RCW103 3.3 


1-3 


11, 12, 14 


RXJ2020. 2+4026 


4 x 10~ 14 


7 x 10 6 


< 5 


< 0.35 


G078.2+2.1 


7-Cyg 1.5 


10 


10, 16 


- cases with no SNR - 
















Geminga 


4 x 10" 13 


3 x 10 4 


< 0.1 


< 1.0 




0.16 


350 


8 


RXJ1856. 5-3754 


1 x lO" 11 




< 4 






0.12 




5,6 


MS 0317-6647 


4 x 10~ 13 




? 


? 
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Table 1. Candidate isolated neutron stars and their properties. X-ray fluxes are from ROSAT and are used to derive a spin-down power 
E except where an alternative estimate is listed. These cases are discussed in the text. Radio flux limits are as given in the references 
and arc a mixture of pulse and point source limits as available. Distances and ages are for the supernova remnants, where one has been 
associated with the neutron star candidate. References: 1. Slane et al. 1997; 2. Brazier et al. 1997; 3. Helfand et al. 1995; 4. Petre et al. 
1996; 5. Walter ct al. 1996; 6. Neuhauser et al. 1997; 7. Frail & Moffctt 1993; 8. Sciradakis 1992; 9. Hailcy & Craig 1995; 10. Brazier et 
al. 1996; 11. Tuohy & Garmirc 1980; 12. Gotthelf et al. 1997; 13. Stockc ct al. 1995; 14. Kaspi et al. 1996; 15. Mereghetti et al. 1996; 16. 
Lorimcr ct al. 1997b. 



listed at this position, suggesting that RXJ0002+6246 could 
be a radio-quiet neutron star associated with the proposed 
remnant, possibly with a pulsation period of 242 ms (Hailey 
& Craig 1995). Because no quantitative limit to the radio 
flux is given, RXJ0002+6246 will not be included in this 
paper as a radio-quiet pulsar candidate. 

2.1.2 RXJ0007. 0+7302 (CTA 1) 

This object lies in the north of the CTA 1 supernova rem- 
nant, a mostly circular radio SNR with a 'blow-out' on its 
north-east side (Pineault et al. 1993). In X-rays, the SNR 
is dominated by non-thermal plerionic emission, interpreted 
as synchrotron emission from the relativistic particles gener- 
ated by a fast pulsar (Slane et al. 1997). Slane et al. propose 
that the point source RXJ0007. 0+7302, which is positioned 
at the centre of the plerion, is the pulsar, and this is sup- 
ported by limits on its optical and radio flux (Brazier et al. 
1997). The flux and spectrum of a persistent, unidentified 
7-ray source coincident with RXJ0007. 0+7302 are consis- 
tent with emission from a pulsar at the age and distance of 
CTA 1, leading Brazier et al. (1997) to propose that this ob- 
ject, like Geminga, is a 7-ray loud, radio-quiet pulsar. If we 
accept this proposal, the spectrum peaks in the GeV region, 
making this a 'Vela-like' rather than 'Crab-like' pulsar, with 
the faint flux in soft X-rays typical of pulsars older than a 
few thousand years. 

2.1.3 RXJ0201. 8+6435 (3C58) 

3C58 (G130.7+3.1) is a young SNR, similar to the Crab SNR 
in radio spectrum and morphology, and is thought to have 
resulted from the supernova of AD 1181 (Clark & Stephen- 
son 1977). Confirmation of a compact X-ray source within 
the plerion (Becker et al. 1982) was provided by Helfand et 
al. (1995). It is still unclear whether a north-south ellipticity 
of the X-ray source means that the object is in fact asso- 
ciated with a radio filament discovered by Frail & Moffett 



(1993) or that the apparent extent is due to the attitude re- 
construction problems known to affect some ROSAT data. 
Helfand et al. (1995) determine that the X-ray luminosity 
of the plerionic SNR emission is (4-7) x 10 ergs -1 , which 
would require the neutron star to have a spin-down power of 
E ~ 10 35 - 4 x 10 36 ergs _1 , bracketing the estimated power 
of 1.5 x 10 36 ergs _1 needed to produce a shock at the position 
of the radio filament (Frail & Moffett 1993). This is much 
smaller than than the power of known radio pulsars with 
such small ages. The ROSAT spectrum and the upper limit 
of 50% on pulsations are best modelled as thermal emission 
from the hot polar caps of a neutron star, not non-thermal 
magnetospheric pulses. The properties of the nebula and the 
compact source can be reconciled within the age constraints 
if the pulsar is relatively slow and with a surface magnetic 
field above 10 13 G (Helfand et al. 1995). 

2.1.4 IE 0820-4247 (Puppis A) 

A compact source close to the centre of Puppis A was first 
discovered in Einstein HRI images (Petre et al. 1982), but 
its point-like nature has only recently been confirmed (Pe- 
tre et al. 1996). Stringent optical and radio limits (Petre et 
al. 1996, Kaspi et al. 1996) rule out most types of X-ray 
sources apart from a neutron star or a BL Lac with weak 
radio emission. The X-ray source shows no evidence of vari- 
ability between observations and is less than 20% pulsed, 
implying that, if this is a neutron star, the soft X-ray emis- 
sion is predominantly thermal rather than magnetospheric, 
i.e. more similar to the Vela pulsar than to the Crab. There 
is no visible plerion. 

2.1.5 IE 1207.4-5209 (G296. 5+10.0) 

G296. 5+10.0 has a symmetric 'barrel' morphology defined 
by two clear, elongated radio arcs; its age is estimated to be 
roughly 20,000 years (Seward & Wang 1988). Using Einstein 
observations of the remnant, Kellett et al. (1987) showed 
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that a compact X-ray source, IE 1207.4-5209, close to the 
geometric centre of the remnant might be an isolated neu- 
tron star. Matsui et al. (1988) strengthened this claim with 
upper limits on the optical luminosity, while Mereghetti et 
al. (1996) have published further details of the effort to iden- 
tify this X-ray source. They find no plausible optical coun- 
terpart down to mv ~ 25 and place a 0.1 mjy limit at 4.8 
GHz. Kaspi et al. (1996) place a tighter 1 mjy pulsation 
limit at 400 MHz. Recent analysis of the ASCA/ROSAT 
spectrum of IE 1207.4-5209 shows that it is thermal and 
can be interpreted in terms of a cooling neutron star (Va- 
sisht et al. 1997). Less than 24% of the ASCA flux is pulsed. 

2.1.6 IE 161348-5055 (RCW 103) 

Einstein HRI observations first revealed IE 161348-5055 
as an unresolved X-ray source within RCW 103 (Tuohy & 
Garmire 1980). With the object lying very close to the cen- 
tre of the SNR in a minimum of the diffuse shell emission, 
an association between the two seemed likely, but searches 
for optical or radio counterparts have not been successful. 
Kaspi et al. (1996) provide 400 MHz and 1.5 GHz pulsed flux 
density limits of 1.5 mjy and 0.1 mjy respectively. Recently, 
the point source has been observed with ASCA and sepa- 
rated spectrally from the bright but softer SNR shell emis- 
sion (Gotthelf et al. 1997). The point source flux is probably 
constant on long time scales, within the errors presented by 
matching an ill-defined spectrum across several instruments. 
The age of the remnant is estimated to be 1,000-3,000 years, 
and the supernova responsible for RCW 103 might have been 
a guest star reported in 134 B.C. (Wang et al. 1986). 

2.1.7 RXJ2020. 2+4026 (G078. 2+2.1) 

The 7-ray source 2EG J2020+4026 (2CG078) has been 
linked with the G078. 2+2.1 SNR by several authors, most 
recently Sturner & Dermer (1995) and Esposito et al. (1996). 
During a detailed study of the 7-ray source, Brazier et 
al. (1996) noted the existence of a single unresolved X-ray 
source, RXJ2020. 2+4026, at the centre of the remnant. The 
X-ray flux is steady and no likely optical or radio counter- 
parts were found in subsequent searches, leaving the possi- 
bility that the X-ray and 7-ray fluxes were from a pulsar in 
the G078. 2+2.1 SNR. Like RXJ0007.0+7302, the proposed 
pulsar is Vela-like, with a small X-ray flux relative to the 
7-ray flux. 

2.2 Isolated candidates 

2.2.1 Geminga 

The discovery of 4.2 Hz pulsations in the enigmatic Geminga, 
first in X-rays (Halpern & Holt 1992) and then in 7-rays 
(Bertsch et al. 1992), confirmed this as a pulsar, albeit one 
with a radio luminosity very much lower than any known 
radio pulsar. Recent reports that Geminga has finally been 
seen as a radio pulsar at very low frequencies (Kuzmin & 
Losovsky 1997, Malofeev & Malov 1997) do not affect the 
limits on its luminosity at higher frequencies, which we use 
here for comparison with the population of radio pulsars. 

The lack of absorption of Geminga's X-ray flux and 
its high 7-ray flux mean that Geminga must be nearby, 
within the approximate range 100 < d < 400 pc (Bertsch 



et al. 1992, Halpern & Ruderman 1993). A recent, marginal 
detection of parallax in Hubble Space Telescope observa- 
tions gives a distance of 120-220 pc (Caraveo et al. 1996), 
in good agreement with the earlier estimates. The best par- 
allax value of 160 pc will be adopted in this paper 

Knowledge of Geminga's spin parameters and distance 
enables a direct comparison with radio pulsars. In particu- 
lar, the spin-down power and age can be derived from the 
rotation period P and its first derivative P in exactly the 
same way as for other pulsars. 

2.2.2 MS 0317-6647 

Among the list of candidate neutron stars, MS 0317-6647 is 
the only one to show signs of long-term variability. It is an 
unusual source, with no optical counterpart and a hard, fea- 
tureless X-ray spectrum that is not well described by simple 
models (Petre et al. 1994) . While it may well be a luminous 
X-ray binary in the spiral galaxy NGC 1313, Stocke et al. 
(1995) discuss the possibility that it is an old neutron star 
in our own galaxy. The variability argues against a rotation- 
powered or cooling neutron star, leaving the option that this 
is an object accreting from the interstellar medium (ISM). 
No specific radio flux limits are available, and this object is 
not included further in this paper. 

2.2.3 RXJ1856. 5-3754 

Walter et al. (1996) announced the discovery that this X- 
ray source was an old neutron star at a distance of about 
100 pc. The source, first detected in the Einstein slew sur- 
vey, is very bright, yet no counterpart was found, and the 
source flux did not vary. The proposed distance to the ob- 
ject was determined from the X-ray spectrum of the object: 
since the source is in the direction of a molecular cloud at 
~ 120 pc, the low Nh measured from the X-ray spectrum 
implied that the source lay in front of the cloud. Walter et 
al. concluded that the object was a nearby, old neutron star, 
perhaps powering its X-ray emission through accretion from 
the interstellar medium. 

Reassessments of RXJ1856. 5-3754 have been published 
by Neuhauser et al. (1997) and Campana et al. (1997). Both 
groups find that their more conservative error circles for the 
X-ray position still exclude all bright optical sources, and 
they agree with the previous conclusion that the object is 
probably a neutron star. However, it has not been possible 
to distinguish between a "middle-aged" pulsar like Geminga 
and an old neutron star accreting from the ISM. Using the 
revised source position given by Neuhauser et al., Walter & 
Matthews (1997) have now identified a faint optical counter- 
part at magnitude 25.6, consistent with a neutron star. The 
distance to the object may be larger than claimed by Walter 
ct al. (1996). The X-ray spectrum is well modelled by black 
body radiation from the surface of a neutron star at 100-170 
pc. However, the X-ray spectra of radio pulsars do not give 
reliable distances under the same assumptions. Observations 
of the parallax and/or proper motion of RXJ1856. 5-3754 
will help to resolve its nature and distance. For the rest of 
this paper we will assume that it is at a distance of 120 pc. 

From the ten objects listed in Table 1 we conservatively 
accept the eight with quantitative radio limits as good neu- 
tron star candidates, six in SNRs and two isolated. The 
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Compact source SNR SNR radius Offset fi v 







(arcmin) 


(arcmin) 




(km s 


RXJ0007.0+7302 


CTA 1 


45 


15 


0.33 


730 


RXJ0201. 8+6435 


3C58 


3 


? 


? 


? 


IE 0820-4247 


Pup A 


25 


6 


0.25 


930 


IE 1207.4-5209 


G296.5+10.0 


33 


3 


0.10 


200 


IE 161348-5055 


RCW103 


10 


0.5 


0.05 


220 


RXJ2020.2+4026 


G078. 2+2.1 


60 


1.5 


0.03 


60 



Table 2. Offsets of pulsar candidates from SNR centres, expressed as fractions 
(3 of the SNR radius. The implied velocities, using the distance and age estimates 
in Table 1, are also given. Supernova sizes and centres are from Green's SNR 
catalogue (1996) and from Petre et al. (1996). The size of the 3C58 remnant is 
for the plerion alone. 



strong case for physical association between the first group 
of X-ray sources and the supernova remnants is seen more 
clearly if we consider the distances from the point sources 
to the SNR centres and the transverse velocities that these 
distances imply. With the exclusion of 3C58, which has no 
known SNR shell, these are listed in Table 2. The velocities 
are entirely consistent with the median of 460 kms -1 given 
by this method for radio pulsars in SNRs (Frail et al. 1994) 
and the median of 300 kms -1 seen in the general popula- 
tion of radio pulsars (Lorimer et al. 1997a). Also listed in 
Table 2 is the angular distance /3 between the pulsar and the 
SNR centre, expressed as a fraction of the SNR radius. The 
values of f3 here are smaller than in most suggested associa- 
tions between radio pulsars and supernova remnants (Kaspi 
1996 & references therein). It is improbable that so many 
unusual X-ray objects would be found near the centres of 
young SNRs unless there is an association. In addition, the 
plerions in CTA 1 and 3C58 point very strongly to active 
sources of relativistic particles in these SNRs. The lack of 
a plerion in the other SNRs, however, does not imply the 
reverse (Bhattacharya 1990). 

For the six candidates in supernova remnants, there is 
therefore strong evidence that they are neutron stars. Of 
the two isolated candidates, Geminga is confirmed as a pul- 
sar and RX J1856. 5-3754 has properties entirely consistent 
with the proposition that it too is a neutron star. For the 
remainder of this paper, it will be assumed that all eight 
candidates are neutron stars. 

3 COMPARING RADIO PULSARS WITH 
X-RAY-SELECTED NEUTRON STAR 
CANDIDATES 

None of the eight accepted objects has been detected as a 
radio source. Could any of them be a normal radio pulsar? 
In Fig. 1 the radio luminosities, assuming 1 steradian beam- 
ing, of radio pulsars listed in the Princeton pulsar catalogue 
(Taylor et al. 1993) are plotted against spin-down power. 
Larger symbols indicate pulsars at distances of less than 3.5 
kpc, directly comparable with the X-ray pulsar candidates. 
The axes demonstrate the rise towards a higher median lu- 
minosity, ~ 300mJykpc 2 at 400 MHz, ~ lOOmJykpc 2 at 
1400 MHz, for energetic pulsars (e.g. Lyne et al. 1985, Tay- 
lor & Stinebring 1986, Tauris & Manchester 1997) and the 
sensitivity limits of radio searches. Using Table 1, the upper 



limits for the seven pulsar candidates and Geminga have also 
been added. Each of the candidate pulsars is shown with an 
order of magnitude uncertainty in the spin-down power: this 
is intended to be illustrative and does not indicate a formal 
uncertainty. Note that for the candidates, both axes scale 
with the square of distance. 

An alternative way to view these data is in terms of pul- 
sar/candidate age, as shown in Figure 2. In this figure, the 
ages for the candidate neutron stars are taken to be those 
of the host SNR, again giving a dependence on distance. 
Increasing the distances in both figures, for example, would 
tend to make the radio limits weaker and the candidate pul- 
sars older but also more powerful. Even taking reasonable 
degrees of uncertainty into account, very large increases in 
distance would be required to bring the radio luminosities of 
the candidate pulsars into the realm of known radio pulsars. 
There are four possible explanations for their non-detection 
as radio sources: 

1. They have intrinsically small radio luminosities or very 
steep radio spectra 

2. The radio beams are directed away from Earth 

3. They have unusual spin parameters untested by current 
pulsar searches. 

4. They are not pulsars. 

Although we have accepted all eight candidates as neu- 
tron stars, it will of course remain possible that one or more 
of the candidates is not a rotation-powered pulsar until each 
of them has been identified conclusively. Weak accretion has 
been suggested as a possible power source in RX J 1856. 5- 
3754 and IE 1207.4-5209. The accretion models fall into 
two types: old neutron stars accreting from the interstellar 
medium and 'anomalous X-ray pulsars' (van Paradijs et al. 
1995) with periods of a few seconds and large X-ray lumi- 
nosities relative to the power available from their spin down. 
The latter objects are thought to be isolated neutron stars 
powering their X-radiation by accretion and to result from 
a high mass X-ray binary (HMXB) which underwent a com- 
mon envelope phase (van Paradijs et al. 1995, Ghosh et al. 
1997). The distances and ages of the four examples known 
give a total of only 20-50 such pulsars in the Galaxy, a 
very low number appropriate for their exotic history. There- 
fore, although these objects are difficult to distinguish from 
other pulsars without the detection of pulsations and are 
radio-quiet, it is unlikely that they explain all the candi- 
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Figure 1. The luminosities of radio pulsars at 400 MHz (top) and 1.4 GHz (bottom) against spin-down power. Larger symbols indicate 
pulsars at distances of less than 3.5 kpc. Restrictive upper limits are included for Geminga and the X— ray-selected candidate neutron 
stars selected from Tabic 1. Dashed horizontal lines illustrate a luminosity of lmjykpc 2 . 



date neutron stars in Table 1. Old neutron stars producing 
detectable X-ray fluxes by accretion from the interstellar 
medium must also be rare because deep searches have failed 
to find any (Manning et al. 1996, Danner 1997a, 1997b). 
This mechanism is a possible explanation for only one of 



the candidate neutron stars, RXJ1856. 5-3754. 

Empirically, the underlying luminosity function of nearby 
radio pulsars flattens below 20 mjy kpc 2 (Lyne et al. 1997), 
giving no evidence for substantial numbers of low luminos- 
ity pulsars. In addition, proposing that the candidate pulsars 
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Figure 2. The luminosities of radio pulsars and candidate neutron stars at 400 MHz (top) and 1.4 GHz (bottom) against pulsar or SNR 
age, where pulsar ages are calculated from r = p/2p and SNR ages are listed in Table 1. Larger symbols show the positions of pulsars 
with distances of less than 3.5 kpc. 



have low radio luminosities or steep spectra leads to a con- 
flict with such population constraints, unless there is some 
mechanism that can suppress the radio emission from young 
pulsars but permits high energy radiation and the formation 
of plerions. This is a particular problem for the candidates in 



young SNRs, since we would expect large numbers of similar 
objects at greater ages. 

It has previously been proposed that pulsars with ex- 
treme parameters, such as unusually long periods or high 
magnetic fields, might have radio luminosities low enough 
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to escape detection. For example, slow (P > 0.3 s) pulsars 
with low magnetic fields might be 'injected' into the popu- 
lation, forming a sub-population of low luminosity pulsars 
(Narayan 1987). However, there is no compelling observa- 
tional evidence for injected pulsars. A number of authors 
have now shown that the results of radio pulsar surveys and 
the small number of pulsar/SNR associations can be ex- 
plained at least as well by a simple, single population as 
by multi-population models (Bhattacharya et al. 1992, Frail 
& Moffett 1993, Lorimer et al. 1993, Gaensler & Johnston 
1995, Lorimer et al. 1997b). 

The simplest explanation for the candidate radio-quiet 
pulsars is that they are radio pulsars whose beams do not 
sweep past the Earth. This requires no new population of 
neutron stars and can accommodate all of the candidates 
and Geminga. However, it has a number of implications for 
pulsar beaming and pulsar statistics, discussed below. 

3.1 Beaming 

Pulsars are complex X-ray sources. The very young, Crab- 
like radio pulsars are visible in soft X-rays as non-thermal, 
highly pulsed, compact sources surrounded by a bright neb- 
ula of synchrotron X-rays. Older pulsars, up to a few xlO 5 
years old, are generally much weaker point sources and usu- 
ally lack a synchrotron nebula. In these pulsars, the non- 
thermal pulses become dominant only above a few keV, 
while the X-rays in the ROSAT (0.1 - 2.4 keV) band are 
thermal radiation from the 10 6 K neutron star surface, repre- 
senting the cooling of the hot stellar interior and bombarding 
of the polar caps by fast particles from the magnetosphere. 
The polar caps of the neutron star are hotter than the rest 
of the surface and are detected as pulses, but because the 
radiation is bent gravitationally, the pulses are spread out 
to give a low degree of modulation even when radiation is 
entirely from the polar caps (Zavlin et al. 1995, Page 1995, 
Yancopoulos et al. 1994). Typically less than 20% of the keV 
flux in detected pulsars forms the characteristically broad, 
smooth pulse, in agreement with the predictions. We will 
therefore assume that this thermal radiation is spread by 
gravitational bending to be visible from all directions. This 
means that our ability to detect a neutron star from its soft 
X-ray emission is unaffected by beaming. 

Thermal emission appears to be responsible for the X- 
rays from the candidate neutron stars, given their low lu- 
minosities, soft spectra and the lack of pulsations below 2 
keV. Where the spectra have been modelled, they have been 
described as blackbody, although not all of them are well 
constrained. 

The beaming fraction of radio pulsars is not well known. 
It is generally accepted that the average beaming frac- 
tion for the whole population is around 20% (e.g. Lyne & 
Manchester 1988). Tauris & Manchester (1997) find that 
it is only 10% but anti-correlates with age, giving a much 
higher beaming fraction for young pulsars. Some authors 
(e.g. Narayan 1987) find a beaming fraction approaching 
100% for young pulsars, boosting their case for injection 
of long-period pulsars into the population. Frail & Moffett 
(1993) derived a value of 61 ± 13% from a deep search for 
(young) radio pulsars in plerions. 

We can provide an independent measure of the beaming 
fraction by considering the pulsars which have been detected 



in X-rays but are radio-quiet. Taking just the pulsar candi- 
dates with quantitative radio flux limits, the six in SNRs 
have ages less than 20,000 yr and distances smaller than 3.5 
kpc. Six known radio pulsars also satisfy these criteria: PSRs 
B0531+21, B0833-45, B1706-44 and B1046-58, all of which 
have been detected at high energies, and PSRs B1737-30 
and B1853+01, which have not. These latter two pulsars 
are relatively low down on the E/d 2 ranking, probably ex- 
plaining why they have not been detected. The number of 
radio-quiet pulsars relative to the total yields a radio beam- 
ing fraction of ~ 50% for young pulsars, assuming that the 
sample is largely complete. We argue below that this is the 
case. The radio-quiet pulsars are clearly inconsistent with 
beaming fractions close to 100%. 

3.1.1 Geminga-like 'y-ray pulsars 

The presence of pulsed 7-rays from a radio-quiet pulsar 
allows us to constrain the emission geometry of the radio 
and high energy beams. At present 7-ray pulses have been 
identified only in Geminga, although the candidate pulsars 
RXJ0007.0+7302 and RXJ2020. 2+4026 are also coincident 
with 7-ray sources. The beaming explanation for radio-quiet 
neutron stars therefore demands a model in which it is pos- 
sible to see the hard X-/7~ray pulses without intersecting 
the radio beam. 

Current data on the high energy emission from pulsars 
is limited by the sensitivity of available instruments. In 7— 
rays, just six radio pulsars have been detected as pulsed 
sources (Thompson et al. 1994, Ramanamurthy et al. 1995, 
Carraminana et al. 1995). The range of pulse shapes, from 
a single broad hump to two widely separated, sharp peaks 
connected by a saddle, can be explained in terms of dif- 
ferent lines of sight across a single, edge-brightened beam 
(Daugherty &: Harding 1996, Romani & Yadigaroglu 1995). 
The wide pulses imply a broad beam unless there are only 
small offsets between the observer, the beam axis and the 
neutron star spin axis. Wider beams are generally prefer- 
able because they do not require such a specific geometry 
and can explain why all of the radio pulsars with highest 
E/d 2 have been detected. 

In order to explain radio-quiet pulsars by beaming, the 
high energy and radio beams cannot be generated in the 
same location in the neutron star magnetosphere. One possi- 
bility (e.g. Romani & Yadigaroglu 1995) is that the high en- 
ergy beams are produced near the light cylinder (co-rotation 
radius) while the radio is beamed along the magnetic axis. 
Efficient production of high energy photons in this model 
depends on a high inclination of the magnetic axis (Romani 
1996) and a pulsar with low magnetic inclination or viewed 
from near the spin axis would not be visible at high energies. 
Approximately two-thirds of the pulsars visible as EGRET 
7-ray sources would be radio-quiet in this model (Yadi- 
garoglu & Romani 1995), although a larger radio beam- 
ing fraction and allowance for deep, targeted radio pulsar 
searches will decrease this figure. Several groups are working 
on unidentified 7-ray sources to look for radio-quiet pulsars. 

Uncertainty over the inclinations measured from ra- 
dio polarisation (Lyne & Manchester 1988, Rankin 1990, 
Manchester 1996) fuels an on-going debate between the 
above model and an alternative set of models with exactly 
the opposite geometric requirements (Daugherty & Harding 
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1996, Sturner & Dermer 1995). These 'polar cap' models 
have been successful in explaining the spectra and 7~ray lu- 
minosities of radio pulsars, but they require that the pulsar 
magnetic and spin axes are approximately aligned in order 
for the emission from a single magnetic pole to produce the 
broad observed pulses. The most recent simulations (Daugh- 
erty & Harding 1996) have more relaxed geometries than 
earlier models. If polar cap models are to explain radio-quiet 
pulsars by beaming, then the radio beams come either from 
a separate region of the magnetosphere or are internal to 
the hollow 7-ray cone. The radio pulses of 7-ray pulsars are 
(so far) always outside the 7-ray pulse, inconsistent with an 
internal radio beam (Daugherty & Harding 1996). In both 
outer gap and polar cap models, therefore, the radio and 
high energy beams of radio-quiet pulsars must be generated 
in different parts of the magnetosphere. 

3.2 The pulsar birth rate 

Pulsar birth rates are usually derived from observations and 
models of the radio population alone. In this section we use 
the X-ray observations of young neutron stars, usually ne- 
glected in such calculations, to constrain the local pulsar 
birth rate independent of radio beaming and luminosity laws. 

As we listed in Section 3.1, 10 neutron stars with ages 
less than 20,000 years and distances below 3.5 kpc have 
been detected in X-rays. This implies a birth rate of 13 
Myr -1 kpc -2 , if pulsars are born close to the Galactic disk 
and the X-rays are not beamed. In order to extrapolate this 
to the whole Galaxy, we assume that the radial distribution 
of pulsars is a Gaussian with a radial scale length of 5 kpc 
(e.g. Lorimer et al. 1993), which gives a Galactic neutron 
star birth rate of 1 every 110 years. Adding two radio pul- 
sars that have not yet been detected in X-rays raises this 
figure to 1 every 91 years. Further allowance for incomplete- 
ness in the X-ray detections can only increase the birth rate 
further. 

The frequency of supernova explosions and the birthrates 
of SNRs and pulsars have been the subject of discussion 
since the late 1960's. Recent calculations give a total (Type 

I plus Type II) Galactic supernova rate of 1 every 40 years 
(Tammann et al. 1994), similar to the rate of one Type 

II every 50-170 years derived from extragalactic supernova 
searches (Cappellaro et al. 1997). A recent estimate for the 
birth rate of radio pulsars is 1 every 60 - 330 years (Lyne 
et al. 1997). Gaensler & Johnston (1995) found that a birth 
rate of one every 85 years gave an excellent match between 
observed and modelled SNR/pulsar associations. 

The birth rate cannot be a factor of 2 higher than we 
have derived from the X-ray neutron star population, or it 
will be in conflict with the (independently derived) super- 
nova rate. It is also unlikely to be much lower than our es- 
timate because our sample is not complete. The birth rates 
for neutron stars, Type II supernovae and radio pulsars are 
therefore similar. We conclude that (a) probably all young 
neutron stars are radio pulsars, (b) more young pulsars are 
visible as X-ray sources than as radio pulsars, and (c) most 
of the young, nearby pulsars have already been discovered. 

Lastly, our result supports the conclusion of Gaensler 
& Johnston (1995) that the small number of known pul- 
sar/SNR associations is a consequence of pulsar beaming 
and luminosity, not a dearth of radio pulsars. 



4 CONCLUSIONS 

This paper has clarified the evidence for radio-quiet pul- 
sars and the implications of such objects. We have listed six 
clear, unresolved X-ray sources in supernova remnants, with 
quantitative radio flux limits and high X-ray to optical flux 
ratios that rule out nearly all types of X-ray source. Most of 
them are the only X-ray point source within their SNR and 
they are all close to their SNR centres, making a strong case 
that they are stellar remnants associated with the SNRs. 
Their transverse velocities are consistent with the velocities 
observed in radio pulsars (Frail et al. 1994, Lorimer et al. 
1997a). We find that it is simplest to explain all of these 
objects, plus two further objects without SNRs, as neutron 
stars. 

The candidate neutron stars have lower radio fluxes 
than would be expected from known radio pulsars of equiv- 
alent age or spin-down power. Reasons for this might in- 
clude extreme spin parameters (e.g because of large mag- 
netic fields) or truly low radio luminosities. However, these 
are not necessary to explain the sources or justified by other 
empirical evidence. The low radio luminosities are most sim- 
ply accommodated in a geometric explanation, in which the 
radio emission is not favourably beamed whereas the soft 
X-rays are dominated by thermal emission from the neu- 
tron star surface and are visible from all directions. The 
relative numbers of radio pulsars and X-ray pulsar candi- 
dates in SNRs gives a crude estimate of ~ 50% for the radio 
beaming fraction. 

Above ~ 2 keV, non-thermal emission from the magne- 
tosphere becomes dominant in known pulsars. The presence 
of high energy radiation without radio pulses implies differ- 
ent emission sites for the two ends of the spectrum. The can- 
didate pulsars in CTA 1 and G078. 2+2.1 coincide with 7-ray 
sources and searches for pulsations in the high energy fluxes 
should be pursued. IE 1207.4-5209, the candidate neutron 
star in G296. 5+10.0, has not shown any evidence for mag- 
netospheric emission. This could mean that the object is a 
cooling neutron star with only weak magnetospheric activ- 
ity (Vasisht et al. 1997), but it could also be a pulsar in 
which both the radio and high energy beams are directed 
away from the Earth. The object's rotation frequency may 
still be discovered from low-level modulations of the thermal 
soft X-rays. 

We have also used the assumption of quasi-isotropic X— 
ray emission to estimate the neutron star birth rate, which 
we find to be at least 13 Myr -1 kpc -2 in the neighbourhood 
of the Sun. The total galactic birth rate is therefore at least 
1 neutron star every ~ 90 years, close to the derived rate of 
Type II supernovae (Cappellaro et al. 1997). We conclude 
that neutron stars are a frequent outcome of supernovae, 
that probably all neutron stars are born as radio pulsar and 
that most young, nearby pulsars have already been discov- 
ered. This is further support for our result that radio-quiet 
pulsars are best explained as unfavourably beamed radio 
pulsars. 
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